Electronic States of Graphene Nanoribbons 
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We study the electronic states of narrow graphene ribbons ("nanoribbons") with zigzag and 
armchair edges. The finite width of these systems breaks the spectrum into an infinite set of bands, 
which we demonstrate can be quantitatively understood using the Dirac equation with appropriate 
boundary conditions. For the zigzag nanoribbon we demonstrate that the boundary condition 
allows a particle- and a hole-like band with evanescent wavefunctions confined to the surfaces, 
which continuously turn into the well-known zero energy surface states as the width gets large. For 
armchair edges, we show that the boundary condition leads to admixing of valley states, and the 
band structure is metallic when the width of the sample in lattice constant units is divisible by 3, and 
insulating otherwise. A comparison of the wavefunctions and energies from tight- binding calculations 
and solutions of the Dirac equations yields quantitative agreement for all but the narrowest ribbons. 

PACS numbers: 73.22-f,73.20-r,73.23-b 



I. INTRODUCTION 

Improvements in the processing of graphite have made 
possible the isolation of two dimensional carbon sheets 
known as graphene I]- The experimental observation 
of the quantum Hall effect in ribbons of this material, 
with widths in the micron Q or submicron [j| range, in- 
dicates unambiguously the two dimensional character of 
the system. The possibility of gating and further pro- 
cessing such graphene sheets into multi-terminal devices 
has opened a new field of carbon-based nanoelectronics, 
where graphene nanoribbons could be used as connec- 
tions in nanodevices. 

The electronic properties of a nanometer scale carbon 
system depends strongly on its size and geometry 0, |j| . 
This is well known in the case of nanotubes, which are 
graphene sheets rolled into cylinders 6]. The geome- 
try dependence is strongly influenced by the bipartite 
character of the graphene lattice. For carbon nanotubes 
the wrapping direction imposes different boundary con- 
ditions on the wavefunction in the different sublattices, 
which determines whether the system is semiconducting 
or metallic. 

In this work we study the electronic states of graphene 
ribbons with different atomic terminations (FigQJ. Us- 
ing tight-binding calculations, we show that the elec- 
tronic properties depend strongly on the size and geom- 
etry of the graphene nanoribbons. We demonstrate that 
the electronic energies and states may be understood in 
terms of eigenvalues and eigenvectors of the Dirac Hamil- 
tonian, which describes the physics of the electrons near 
the Fermi energy of the undoped material. 

We now summarize our results. We find that for 
nanoribbons with zigzag edges, the correct boundary con- 
dition is for the wavefunction to vanish on a single sub- 
lattice at each edge. In this case the nanoribbon has 
confined electronic states with wavefunctions that involve 
sites on both sublattices, and is extended across the sys- 
tem. In addition, there are surface states strongly lo- 
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FIG. 1: (Color online)The lattice structure of a graphene 
sheet. The primitive lattice vectors are denoted by a and b. 
Top and bottom are zigzag edges, left and right are armchair 
edges. Atoms enclosed in the vertical (horizontal) rectangle 
represent the unit cell used in the calculation of nanoribbons 
with zigzag (armchair) edges. The length of the nanoribbons, 
L, as function of the number of atoms, N, in the unit cell is 
also indicated. 



calized near the edges which are non-vanishing only on 
a single sublattice. For armchair edges, the appropriate 
boundary condition is for the wavefunction to vanish on 
both sublattices at the edges. This can be achieved by 
admixing states from both Dirac points. In this case we 
find that the electronic structure depends critically on 
the nanoribbon width, with the system being metallic 
for nanoribbons of width L = 3Afao, where M an inte- 
ger and do the graphene lattice constant, and insulating 
otherwise. 
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II. MODEL HAMILTONIAN 

In graphene the carbon atoms crystallize in a honey- 
comb structure whose primitive lattice vectors are a = 
ao(l,0) and b = a (l/2, \/3/2). The lattice is bipar- 
tite and there are two atoms per unit cell, denoted by 
A and B, located at (0,0) and at d = ao(0, l/\/3). In 
the simplest model, the carriers move in the x-y plane 
by hopping between the p z orbitals of the carbon atoms. 
A tight-binding model with only nearest neighbor hop- 
ping t leads to a Hamiltonian with Dirac points at the 
six corners of the Brillouin zone, only two of which are 
inequivalent. We take these to be K = ancl 

= a^(~5' 71 )■ Wavefunctions can be expressed via 
the k • P approximation 0, @ m terms of envelope func- 
tions [^A(r), V , s( r )] an d Wa( v )^'b( v )\ f° r states near 
the K and K' points, respectively, which may be com- 
bined into a 4- vector ^> = (tpA, i>B, ~ — ^b) 0- This 
satisfies a Dirac equation H^> = e^, with 
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where 7 = ^/H>t/2. Note that k denotes the separation 
in reciprocal space of the wavefunction from the K (K') 
point in the upper left (lower right) block of the Hamil- 
tonian. 

The bulk solutions of Hamiltonian are well- 
known 0- The eigenstates retain their valley in- 
dex as a good quantum number and the wavefunc- 
tions, with energies e = ±700 |k|, may be written 
as [e lWr e~ ie ^ 2 ,Te ikr e l6k / 2 ,0,0} for the K valley , and 
[0, 0, e lkr e l9k / 2 , ± e lWr e- ie ^ 2 } for the K' valley. Here (9 k = 
arctan k x /k v . Note that a solution to the Dirac equa- 
tion (ipA, ipB, — — w hh energy e has a particle- 
hole conjugate partner [ljj {ipAi ~ tpBi — ^'ai ^'b) with 
energy — z. Because of this, the eigenstates of E 
must be normalized on each sublattice separately 
/dr[|W(r)| 2 + |^(r)| 2 ] = 1/2, for M = A,B. 
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FIG. 2: (Color online) Examples of energy bands for a 
graphene nanoribbon with periodic boundary conditions in 
one direction, k is the wavevector parallel to the nanoribbon 
edge, measured with respect the center of the Brillouin cen- 
ter, (a) Ribbon terminated in zigzag edges with 56 atoms in 
the unit cell. The dispersionless states correspond to confined 
surface states. The band structures of insulating and metallic 
armchair nanoribbons are plotted in (b) and (c) respectively. 



III. ZIGZAG NANORIBBONS 

The geometry of a nanoribbon with zigzag edges is 
illustrated on the top and bottom edges of Fig. ^ It is 
interesting to note that the atoms at each edge are of the 
same sublattice (A on the top edge of Fig.^and B on the 
bottom edge). In Fig. ^ we also show the unit cell used 
in the tight-binding calculations of the zigzag ribbons, 
containing N/2 A- type atoms that alternate along the 
unit cell with N/2 B-type atoms. The total width of 
the nanoribbon is L = -jV3ao- We impose periodic 
boundary conditions along the direction parallel to the 
edge. In our discussions we will assume that the edges 
lie along the y direction, so in in the discussion of the 
zigzag nanoribbons, the coordinate axes in Fig.^wiU be 



rotated by 90°, and the eigenstates are proportional to 
e lky . In Fig. [21 we plot an example of the band structure 
of a nanoribbon with zigzag edges. The finite width of 
the ribbon produces confinement of the electronic states 
near the Dirac points. In Fig.[3]we plot the energy of the 
first three confined states, at k = K v , as a function of 
the nanoribbon width. The two bands of dispersionless 
localized surface states 0,0,0,0 that occur between 
K y and K' y in Fig. |3[a) are also affected by the finite 
width: they admix, and the two bands are slightly offset 
from zero. The dependence of the electronic states on the 
width of the nanoribbon may be understood in terms of 
eigenstates of the Dirac Hamiltonian with appropriate 
boundary conditions: setting the wavefunction to zero 
on the A sublattice on one edge, and on the B sublattice 
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FIG. 3: (Color online) Calculated confined state energies at a 
Dirac point versus the nanoribbon width, in a zigzag nanorib- 
bon. The dots are tight binding results, and the crosses are 
the results of the k ■ P approximation. 



for the other. We can understand the lines of vanishing 
wavefunction to be lattice sites that would lie just beyond 
the edges if bonds had not been cut to form them. 

For the continuum description, we begin by rotating 
the wavevectors in Eq. ^ k x — > k y , k y — > — k x so that 
the zigzag edge lies along the y, and our wavefunctions 
exist in the space < x < L. Translational invariance 
in the y direction guarantees the wavefunctions can be 
written in the form ^(^(r) = e lkyV '<^ M (')(x). For the 
K (K') valley the wavefunctions obey 



(2) 



with e — e/(-fao). It is easy to see if one solves the 
equations for <f>g and <p' A , the remaining wavefunctions 
are determined by 



e<p B 



e<t>' A 



(id x - iky)4> A 
{-id x + ik v )(j>' B . 

The general solutions of Eq. have the form 

<j>p{x) = Ae zx + Be~ zx , 



(3) 



(4) 

with z = \Jk 2 — e 2 , which can be real or imaginary. 

For the zigzag nanoribbon, we meet the boundary con- 
dition for each type of wavefunction separately: 

4> A (x = 0) = cl>' A {x = 0) = ( f> B ^ = L )= ( t>'B(x = L)=0 . (5) 

These conditions leads to a transcendental equation for 
the allowed values of z, 
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ky -)- Z 



(6) 



Eq.© supports solutions with real values of z = k for 
k y > ky = l/L, which correspond to the surface states. 



These have energies ± 



h 2 



k 2 , and are linear combina- 



tions of states localized on the left and right edges of the 
ribbon. For large values of k y , k — ► k y and the surface 
states become decoupled. For k y < there are no states 
with real z that can meet the boundary conditions, so 
surface states are absent. For values of k y in the range 
< k y < ky, the surface states are so strongly admixed 
that, as we show below, they are indistinguishable from 
confined states. 

For pure imaginary z — ik n , the transcendental equa- 
tion becomes 



v tan (k n L) 



(7) 



and for each solution k n there are two confined states 
with energies e — ± A Ik 



k 2 and wavefunctions 



sin(fc n x) 
cos(fc„a;) + k y 



sm(k n x)) 



(8) 



Here the index n indicates the number of nodes of the 
confined wavefunction. Interestingly, for values of k y 
larger than ky, Eq.Q does not support nodeless solu- 
tions, indicating the existence of surface states in this 
region of reciprocal space. The critical value ky is the 
momentum where the lowest energy solution of the tran- 
scendental Eq. El changes from pure real to pure imagi- 
nary, and the energy is equal to ±|fc^|. 

In order to analyze the accuracy of the k • P approxi- 
mation for describing the electronic properties of carbon 
nanoribbons, in Fig. |3| we plot the energies of the three 
lowest confined states of a zigzag nanoribbon as a func- 
tion of its width, both from the tight-binding approach 
and from our solutions to the Dirac equation. It is ap- 
parent that the two approaches match quite well, even 
for rather small widths (~ 35 A). 

In Fig. we plot the squared wavefunction for the low- 
est energy state of a zigzag nanoribbon as obtained in the 
tight binding approach. Fig. 0{a) corresponds to k y = 
(k = — 27r/3ao with respect the center of the Brillouin 
zone), and Fig. g[b) to k y = 0.02 x 27r/3a . The first 
case corresponds to a nodeless confined state, and we 
find the wavefunction is described nearly perfectly by 
Eq.©, whereas the second case is the expected linear 
combination of surface state wavefunctions that decay 
exponentially from the edges as exp (— kx). 



IV. ARMCHAIR NANORIBBONS 

The geometry of a nanoribbon with armchair edges is 
illustrated on the left and right edges of Fig.Q] along with 
the unit cell used in the corresponding tight binding cal- 
culations. In this orientation the width of the nanoribbon 
is related to the number of atoms in the unit cell through 
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FIG. 4: [Color online) Squared wavefunction for the state 
closest to zero energy for a zigzag nanoribbon, as obtained 
from tight binding calculations. The width of the ribbon is 
L = UVSa . (a) k = -27r/3a , and (b) k = (-27r/3o )(l - 
0.02). Both are measured from the center of the Brillouin 



the expression L = a^. Here the edge runs along the 
y direction, and no rotation of the figure is needed to 
represent our calculations. 

The electronic properties of armchair nanoribbons de- 
pend strongly on their width. In Fig.^b) and (c) we plot 
two examples of band structures of armchair nanorib- 
bons. One sees that in the latter figure there is a Dirac 
point, leading to metallic behavior for a non- interacting 
model, whereas the former is a band insulator. In general 
we find that armchair nanoribbons of width L = 3ATao, 
with M integral, are metallic, whereas all the other cases 
are insulators. The energy of the confined states also be- 
have in a discontinuous way with respect to the width 
of the ribbon. In Fig. [5] we plot the energy of the lowest 
(squared) energy confined states at the center of the Bril- 
louin zone as a function of the nanoribbon width. In the 
inset of this figure we see that the separation in energy 
between confined states is also strongly dependent on the 
number of atoms in the unit cell. 

As in the case of the zigzag nanoribbons this behavior 
may be understood in terms of eigenstates of the Dirac 
Hamiltonian with the correct boundary conditions. In 
Fig.n° ne may see that the termination consists of a line 
of A-B dimers, so it is natural to have the wavefunction 
amplitude vanish on both sublattices at x = and x = L. 




FIG. 5: (Color online) Calculated lowest energy confined 
states at the center of the Brillouin zone versus the nanorib- 
bon width, for an armchair nanoribbon. The dots correspond 
to the tight binding results and the crosses are the results of 
the k • P approximation. In the inset we plot the six lowest 
energy confined states for three different widths. The k • P 
results are slightly shifted to the right for clarity. Note that 
the for L = 24ao the k • P results are doubly degenerate. 



To do this we must admix valleys, and require 

M (i = O) = 0^(x = O) 

M* = L) = ^(x = L)e iAKL , 

with AK = With these boundary conditions the 

general solutions of the Dirac equation are planes waves, 



4>b(x) 



Ak n x 



and 4>' B (x) = e ~ iknX . 



(9) 



The wavefunctions on the A sublattice may be obtained 
via Eq. |3| The wavevector k n satisfies the condition 



so that 



e 2ik n L _ e iAK L 



2tt 

2k n L = — j + 2ir n 
3 



with n an integer and j = 0, ±1, determined by 
-=3M + j 



(10) 



(11) 



(12) 



for an integer M. Thus for armchair nanoribbons the 
allowed values of k n are 



i. 4 ( +3 

kn = 7T in ± — 

N \ 3 



(13) 
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FIG. 6: (Color online) Squared wavefunction of the state 
with energy closest to zero for an armchair ribbon of width 
(a) L = 24ao and (b) L — 25ao, as obtained from tight binding 
calculations. 



with energies k^ + ky. Note that this is in contrast 

to the zigzag nanoribbon for which the allowed values of 
k n depend on k y . For a width that is a multiple of 3ao, 
the allowed values of k„, k n = win/N, create doubly 
degenerate states for \n\ > 0, and allow a zero energy 
state when k y — > 0. Nanoribbons of widths that are not 
multiples of three have nondegenerate states and do not 
include a zero energy mode. Thus these nanoribbons are 
band insulators. The quality of the k • P approximation 
for describing the electronic states of armchair nanorib- 
bons is reflected in Fig. [5] where the energies of the con- 
fined states obtained by diagonalizing the tight binding 



Hamiltonian and by solving Eq. l|13|) are compared. The 
quantitative agreement is apparent for all but the nar- 
rowest ribbons, where one does not expect the k • P to 
work well. 

The admixing of different valley states to meet the 
boundary condition means that the wavefunction will os- 
cillate with period 2n/AK 0- This behavior can explic- 
itly be seen in Fig. |SJ which illustrates the squared wave- 
function from the tight binding calculation. The short 
oscillation in the wavefunctions has exactly the period 
expected for the valley mixing we introduced to meet the 
boundary conditions. In the case of Fig.UJa), the ribbon 
is of width is L = 24ao and k n = 0, so that the energy is 
zero and there is no confinement effect on the form of the 
wavefunction. For a ribbon of width L = 25ao [Fig.|BJb)], 
k n is non-zero and one sees a long- wavelength oscillation 
whose period is related to the value of k n . 

V. CONCLUSION 

In this paper we studied eigenstates and eigenenergies 
of graphene nanoribbons using the Dirac equation with 
appropriate boundary conditions, and compared the re- 
sults to those of tight binding calculations. We found 
that except for the narrowest ribbons, the agreement was 
quantitative. Zigzag nanoribbons support surface states 
which go to zero energy in the limit of wide ribbons, 
and can only be found in a k y interval between the Dirac 
points. Armchair nanoribbons have no surface states, but 
in spite of their finite size they have zero energy states for 
appropriately chosen widths, so that the system oscillates 
between insulating and metallic behavior as the width 
changes. Our results show that the continuum descrip- 
tion of graphene may be used in quantitative analysis of 
this system for all but the most narrow systems. 
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